NSF-KITP-08-38 
PUPT- 2261 
arXiv:0803.3295 [hep-th] 



Building an AdS/CFT superconductor 

Sean A. Hartnoll 11 , Christopher P. Herzog" and Gary T. Horowitz' 1 

KITP, University of California 
Santa Barbara, CA 93106, USA 

" Department of Physics, Princeton University 
Princeton, NJ 085 U, USA 

^ Department of Physics, U CSB 
Santa Barbara, CA 93106, USA 

hartnollOkitp . ucsb . edu , cpherzog@princeton . edu , garySphysics . ucsb . edu 



Abstract 

We show that a simple gravitational theory can provide a holographically dual de- 
scription of a superconductor. There is a critical temperature, below which a charged 
condensate forms via a second order phase transition and the (DC) conductivity be- 
comes infinite. The frequency dependent conductivity develops a gap determined by 
the condensate. We find evidence that the condensate consists of pairs of quasiparticles. 



1 Introduction 



A most remarkable result to emerge from string theory is the AdS/CFT correspondence 
[1], which relates string theory on asymptotically anti de Sitter spacetimes to a conformal 
field theory on the boundary. In recent years, it has become clear that this holographic 
correspondence between a gravitational theory and a quantum field theory can be extended 
to describe aspects of nuclear physics such as the results of heavy ion collisions at RHIC 
[2] and to certain condensed matter systems. Phenomena such as the Hall effect [3] and 
Nernst effect [4, 5, 6] have dual gravitational descriptions. One can ask if there is a dual 
gravitational description of superconductivity. 

Conventional superconductors, including many metallic elements (Al, Nb, Pb, ...), are 
well described by BCS theory [7]. However, basic aspects of unconventional superconduc- 
tors, including the pairing mechanism, remain incompletely understood. There are many 
indications that the normal state in these materials is not described by the standard Fermi 
liquid theory [8] . We therefore hope that a tractable theoretical model of a strongly coupled 
system which develops superconductivity will be of interest. Several important unconven- 
tional superconductors, such as the cuprates and organics, are layered and much of the 
physics is 2+1 dimensional. Our model will also be 2+1 dimensional. 

To map a superconductor to a gravity dual, we introduce temperature by adding a 
black hole [9] and a condensate through a charged scalar field. To reproduce the super- 
conductor phase diagram, we require a system that admits black holes with scalar hair at 
low temperature, but no hair at high temperature. While Hertog has shown that neutral 
AdS black holes can have neutral scalar hair only if the theory is unstable [10], Gubser has 
recently suggested that a charged black hole will support charged scalar hair if the charges 
are large enough [11]. We consider a simpler version of Gubser's bulk theory (in which the 
black hole can remain neutral) and show that it indeed provides a dual description of a 
superconductor. 

2 The model: condensing charged operators 

We start with the planar Schwarzschild anti-de Sitter black hole 



ds 1 = -f(r)dt 2 + 



dr 2 



+ r 2 (dx 2 + dy 2 ), 



(1) 



f(r) 



where 




r 2 M 

J2~—- 



(2) 
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L is the AdS radius and M determines the Hawking temperature of the black hole: 

(3) 

This black hole is 3+1 dimensional, and so will be dual to a 2+1 dimensional theory. In 
this background, we now consider a Maxwell field and a charged complex scalar field, with 
Lagrangian density 1 

C = - l -F ah F ab -Vm)-\d^ -iA^\ 2 . (4) 

For simplicity and concreteness, we will focus on the case 

2|^| 2 

vm) = --y-- (5) 

Although the mass squared is negative, it is above the Breitenlohner-Freedman bound [12] 
and hence does not induce an instability It corresponds to a conformally coupled scalar in 
our background (1) and arises in several contexts in which the AdS^/CFT^ correspondence 
is embedded into string theory For instance, the truncation of M theory on AdS± x S 7 
to N = 8 gauged supergravity has scalars and pseudoscalars with this mass, dual to the 
bilinear operators tr<I>( / <3? <7 ) and ti^^ 1 ^- 1 ) in the dual M = 8 Super Yang-Mills theory, 
respectively. However, we should note that our Lagrangian (4) has not been obtained from 
M theory. We expect that our choice of mass is not crucial, and qualitatively similar results 
will hold, e.g., for massless fields. 

We will work in a limit in which the Maxwell field and scalar field do not backreact on 
the metric. This limit is consistent as long as the fields are small in Planck units. (Recall 
that in the analogous case one dimension higher, only Yang- Mills states with energy of order 
N 2 have finite backreaction in the bulk.) Alternatively, this decoupled Abelian-Higgs sector 
can be obtained from the full Einstein-Maxwell-scalar theory considered in [11] through a 
scaling limit in which the product of the charge of the black hole and the charge of the 
scalar field is held fixed while the latter is taken to infinity. Thus we will obtain solutions of 
non-backreacting scalar hair on the black hole. As we shall see, our simple model captures 
the physics of interest. 

Taking a plane symmetric ansatz, ^ = Vl/(r), the scalar field equation of motion is 

„ /f 2\ , f 2 2 



introducing a gauge coupling 1/e 2 in front of the \F\ 2 term in the action is equivalent to rescaling the 
fields * — > and — ► eA^. Setting e = 1 is a choice of units of charge in the dual 2+1 theory. 
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where the scalar potential A t = <£. With A r = A x = A y = 0, the Maxwell equations imply- 
that the phase of *S> must be constant. Without loss of generality we therefore take \E r to 
be real. The equation for the scalar potential $ is the time component of the equation of 
motion for a massive vector field 

2 2* 2 

&> + U' _£_L$ = , (7) 
r } 

where 2^ 2 is the, in our case, r dependent mass. The charged condensate has triggered 
a Higgs mechanism in the bulk theory. At the horizon, r = Vq, for & dt to have finite 
norm, <E> = 0, and (6) then implies ^ = — 3ro i $'/2. Thus, there is a two parameter family 
of solutions which are regular at the horizon. Integrating out to infinity, these solutions 
behave as 

* = — + — + •••• (8) 



and 



$ = M -- + ■■■ ■ (9) 

r 



For both of these falloffs are normalizable [13], so one can impose the boundary condition 
that either one vanishes. 2 After imposing the condition that either or ^( 2 ) vanish we 
have a one parameter family of solutions. 

It follows from (7) that the solution for $ is always monotonic: It starts at zero and 
cannot have a positive maximum or a negative minimum. Note that even though the field 
equations are nonlinear, the overall signs of $ and ^ are not fixed. We will take to be 
positive and hence have a system with positive charge density. The sign of ^ is part of the 
freedom to choose the overall phase of VP. 

Properties of the dual field theory can be read off from the asymptotic behavior of the 
solution. For example, the asymptotic behavior (9) of 3> yields the chemical potential fj, and 
charge density p of the field theory. The condensate of the scalar operator O in the field 
theory dual to the field ^ is given by 

<O i ) = V2* (<) , i = 1,2 (10) 

with the boundary condition = 0. The \/2 normalization simplifies subsequent 

formulae, and corresponds to taking the bulk-boundary coupling \ f d?x ((5^ + O^f). Note 
that 0{ is an operator with dimension i. From this point on we will work in units in which 



2 One might also imagine imposing boundary conditions in which both and are nonzero. How- 
ever, if these boundary conditions respect the AdS symmetries, then the result is a theory in which the 
asymptotic AdS region is unstable [14]. 
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the AdS radius is L = 1. Recall that T has mass dimension one, and p has mass dimension 
two so (Oi)/T l and p/T 2 are dimensionless quantities. 

An exact solution to eqs (6,7) is clearly \F = and = p — p/r. It appears difficult 
to find other analytic solutions to these nonlinear equations. However, it is straightforward 
to solve them numerically. We find that solutions exist with all values of the condensate 
(O). However, as shown in figure 1, in order for the operator to condense, a minimal ratio 
of charge density over temperature squared is required. 




Figure 1: The condensate as a function of temperature for the two operators 0\ and 02- 
The condensate goes to zero at T = T c oc p 1 / 2 . 

The right hand curve in figure 1 is qualitatively similar to that obtained in BCS theory, 
and observed in many materials, where the condensate goes to a constant at zero tempera- 
ture. The left hand curve starts similarly, but at low temperature the condensate appears 
to diverge as T -1 / 6 . However, when the condensate becomes very large, the backreaction 
on the bulk metric can no longer be neglected. At extremely low temperatures, we will 
eventually be outside the region of validity of our approximation. 

By fitting these curves, we see that for small condensate there is a square root behaviour 
that is typical of second order phase transitions. Specifically, for one boundary condition 
we find 

(d) « 9.3T C (1-T/T C ) 1 / 2 , as T T c , (11) 
where the critical temperature is T c « 0.226/? 1 / 2 . For the other boundary condition 

<0 2 > ~ 144T 2 (l-T/r c ) 1 / 2 , as T^T C , (12) 

where 0.118/9 1 / 2 . The continuity of the transition can be checked by computing 

the free energy. Finite temperature continuous symmetry breaking phase transitions are 
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only possible in 2+1 dimensions in the large N limit (i.e. the classical gravity limit of our 
model), where fluctuations are suppressed. These transitions will become crossovers at finite 
N. 

Thus for T < T c a charged scalar operator, (0±) or {O2), has condensed. It is natural 
to expect that this condensate will lead to superconductivity of the current associated with 
this charge. 



3 Maxwell perturbations and the conductivity 

We now compute the conductivity in the dual CFT as a function of frequency. As a first 
step, we need to solve for fluctuations of the vector potential A x in the bulk. The Maxwell 
equation at zero spatial momentum and with a time dependence of the form e~ lu)t gives 



K + ^A' x +^-—)A x = 0. (13) 

To compute causal behavior, we solve this equation with ingoing wave boundary conditions 
at the horizon [16]: A x oc J-W 3r o. The asymptotic behaviour of the Maxwell field at large 
radius is seen to be 

4 (i) 

A x = AM + ^- + ... (14) 

The AdS/CFT dictionary tells us that the dual source and expectation value for the current 
are given by 

A x = 4°) , (J x ) = 4 1 ) . (15) 
Now from Ohm's law we can obtain the conductivity 

a{u) = (JA = = JJdA = . (i6) 

K ' E x Ax uA x Wj4 (o) 1 ; 

In figure 2 we plot the frequency dependent conductivity obtained by solving (13) nu- 
merically. The horizontal line corresponds to temperatures at or above the critical value, 
where there is no condensate. The fact that the conductivity in the normal phase is fre- 
quency independent is a characteristic of theories with AdS& duals [15]. The subsequent 
curves describe successively lower values of the temperature (for fixed charge density). We 
see that as the temperature is lowered, a gap opens. The gap becomes increasingly deep 
until the (real part of the) conductivity is exponentially small. 

There is also a delta function at uo = which appears as soon as T < T c . This can be 
seen by looking at the imaginary part of the conductivity. The Kramers-Kronig relations 
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Figure 2: The formation of a gap in the real, dissipative, part of the conductivity as the 
temperature is lowered below the critical temperature. Results shown for both the 0\ 
operator (left) and the O2 operator (right). There is also a delta function at u = 0. The 
rightmost curve in each plot corresponds to T/T c = .0066 (left) and T/T c = .0026 (right). 

relate the real and imaginary parts of any causal quantity, such as the conductivity, when 
expressed in frequency space. They may be derived from contour integration and the fact 
that functions that vanish when t < 0, when transformed to frequency space only have poles 
in the lower half plane (in our conventions). One of the relations is 

x r , m 1^ f°° Re[<7(a/)]dw' 

From this formula we can see that the real part of the conductivity contains a delta function, 
Re[c(u;)] = tt5(uj), if and only if the imaginary part has a pole, Im[cr(u;)] = 1/lo. One finds 
that there is indeed a pole in Im[<r] at lo = for all T < T c . The superfluid density is the 
coefficient of the delta function of the real part of the conductivity 3 

Re[a(uj)] ~ 7rn s S(uj) . (18) 

By (17), n s is also the coefficient of the pole in the imaginary part Im[cr(u;)] ~ n s /uj as 
co —>■ 0. We find that the superfluid density vanishes linearly with T c — T: 

n s w d (T c — T) as T T c , (19) 

where C\ = 16.5 for the 0\ theory while C2 = 24 for the O2 theory. 



3 The superfluid density is usually defined as the coefficient of S(u>) multiplied by the mass of the electron. 
In simple two fluid models, this density is related to the London magnetic penetration depth, n s — 1/Atv\ 2 l . 
Our scaling (19) thus implies \l ~ (T c — T) -1 / 2 , consistent with Landau-Ginzburg theory. 
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In figure 3 we rescaled the small T/T c plots of figure 2 by plotting the frequency in 
units of the condensate rather than the temperature. The curves tend to a limit in which 
the width of the gap is proportional to the size of the condensate. The differing shapes 
of the plots in figure 3 are precisely what is expected from type II and type I coherence 
factors, respectively [7]. Type II coherence suppresses absorption near the edge of the gap, 
explaining the slower rise of Refer] in the left hand plot. It is possible that this difference is 
due to the operator 0\ being a pair of bosons and O2 a pair of fermions, as in the case of 
AdS A x S 7 . 




< Oi > V < o 2 > 



Figure 3: The gap at small T/T c , with the frequency normalised in terms of the condensate. 
On the right the gap is finite, but since {0\) becomes large at small T, the gap on the left 
is also becoming large. The dashed curve on the left plot is (22). 

The Ferr ell- Glover sum rule states that f Refer] cL; is a constant independent of temper- 
ature. Thus the area missing under the curve Refer] due to the gap must be made up by the 
delta function at cj = 0. That Refer] exceeds the value one in figure 3 (right) implies then 
that the superfluid density n s must be correspondingly reduced for the O2 system compared 
with the 0\ system for T <C T c . 

We can also compute the contribution of the normal, non-superconducting, component 
to the DC conductivity. Let us define 

n n = lim Re[a(u)] . (20) 

From our numerics we obtain 

n n ~ e ~ A / T , for (21) 

where we have A = (0\)/2 and A = \/ (O2) /2. Numerically this factor of 1/2 is accurate 
to at least 4%. From (21), A is immediately interpreted as the energy gap for charged 
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excitations. The gap we found previously in the frequency dependent conductivity was 2A. 
The extra factor of two is expected if the gapped charged quasiparticles are produced in 
pairs, suggesting that there is a 'pairing mechanism' at work in our model. Our results for 
A are suggestive of strong pairing interactions. Note that for figure 1 (right) at T = we 
find 2A ~ 8.4T C , which might be compared with the BCS prediction 2A ~ 3.54T C . The 
larger value is what one expects for deeply bound Cooper pairs. Indeed, in our other model, 
figure 1 (left), we see that A actually diverges at low T. 

Finally in figure 4 we plot the imaginary part of the conductivity in this limiting, low 
temperature, regime. Again we see that the curves rapidly approach a limiting curve. We 
also see the advertised pole at u = 0. 




Figure 4: The imaginary part of the conductivity at small T/T c , with the frequency nor- 
malised in terms of the condensate. The analytic expression (22) is also shown on the left, 
but is indistinguishable from the numerics. 

It is natural to ask if one can reproduce this limiting low temperature behavior by 
just taking M = in our background metric, and considering our matter fields in anti de 
Sitter space. One problem is that there are no solutions to the field equations (6,7) which 
are smooth on the horizon of the Poincare patch. Nevertheless, for the 0\ case, we have 
observed numerically that at low temperatures, ^ ~ {0\)/\/2r. Taking M — > where 
/ « r 2 , (13) can be solved exactly to yield A x = A x ^ exp(±y / (Oi) 2 — oj 2 /r). This exact 
result then produces the nonzero conductivities 

Re[ a ] = V^ 2 ~ (°i) 2 for u> (o 1 ) , lm[a] = VA^iEZZ for w < (d) , (22) 

UJ UJ 

via (16). The curves on the left hand sides of figures 3 and 4 are well approximated by the 
conductivity (22). We have included (22) as a dashed curve in these plots. 
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4 Discussion 



We have shown that a simple 3+1 dimensional bulk theory can reproduce several properties 
of a 2+1 dimensional superconductor. Below a second order superconducting phase transi- 
tion the DC superconductivity becomes infinite and an energy gap for charged excitations 
is formed. 

There are many extensions of this model that we hope to consider elsewhere: 1) By 
probing the system with spatially varying fields and an external magnetic field, one can 
compute the superconducting coherence length and penetration depth, respectively. 2) One 
would like to consider a wider class of models by allowing for more general masses for the 
charged scalar field. 3) One should study the effects of backreaction on the bulk spacetime 
metric. 4) Perhaps the most interesting question is to understand the 'pairing mechanism' 
in field theory that leads to a condensate in these systems. 
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